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Abstract 

Using the Selberg trace formula, we show that for a hyperbohc 
2-orbifold, the spectrum of the Laplacian acting on functions deter- 
mines, and is determined by, the following data: the volume; the total 
length of the mirror boundary; the number of conepoints of each or- 
der, counting a mirror corner as half a conepoint; and the number 
of primitive closed geodesies of each length and orientability class, 
counting a geodesic running along the boundary as half orientation- 
preserving and half orientation-reversing, and discounting imprimi- 
tive geodesies appropriately. This implies that Laplace-isospectral 
hyperbolic 2-orbifolds determine equivalent linear representations of 
Isom(i7^), and are isospectral for any natural operator. 



1 Statement 

We consider compact hyperbolic 2-orbifolds M, not necessarily connected. 
Denote the eigenvalues of the Laplacian acting on functions on M by 

= Ao < Ai < .... 
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We call the sequence (Aq, Ai, . . .) the Laplace spectrum of M. If two spaces 
have the same Laplace spectrum we call them Laplace-isospectral. (Note that 
we don't simply call them 'isospectral', because this term is used in different 
ways by different authors.) 

Our goal here will be to prove: 

Theorem 1. Let M be a compact hyperbolic 2-orbifold, not necessarily con- 
nected. The Laplace spectrum of M determines, and is determined by, the 
following data: 

1. the volume; 

2. the total length of the mirror boundary; 

3. the number of conepoints of each order, counting a mirror corner as 
half a conepoint of the corresponding order; 

4- the number of closed geodesies of each length and orientability class, 
counting a geodesic running along the boundary as half orientation- 
preserving and half orientation-reversing , and counting the k-fold iter- 
ative of a primitive geodesic as worth I. of a primitive geodesic of the 
same length and orientability. 

Of course the Laplace spectrum determines other data as well, for example 
the number of connected components. The data we list here determine those 
other data, since they determine the spectrum. 

Theorem [1] can be recast less picturesquely as follows. Associated to a 
2-orbifold M is a linear representation p^g of Isom(if^) on functions on the 
frame bundle M of M. Associated to this representation is its character Xm^ 
a function on the set of conjugacy classes of Isom(iJ^). The geometrical data 
listed in Theorem [1] are just a way of describing geometrically the information 
conveyed by the character Xm- Selberg tells us that the character determines 
the Laplace spectrum of M. This is a very general fact. What Theorem [1] 
tells us that is special is that for hyperbolic 2-orbifolds, we can get back from 
the Laplace spectrum to the character. 

Once we have the character, we get by general principles the linear equiv- 
alence class of the representation pj^, hence the title, 'Laplace-isospectral 
hyperbolic 2-orbifolds are representation-equivalent'. We also get the spec- 
trum of any natural operator on any natural bundle. We will discuss these 
matters further in IH] below; for now we concentrate on Theorem [TJ 
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2 Plan 



This paper is rather longer than you might think it needs to be, even if 
you disregard the large Appendix. As we observed in [5], this theorem for 
orbifolds is a short step from the result about manifolds proven there. We 
carry on at such length because we are hoping to sell readers on the usefulness 
of the counting kernel technique, and to provide a general background on 
Selberg's methods for those not familiar with them. We also take some 
excursions which we hope will prove interesting. 

We will begin with some examples of Theorem [1] in action; fill in the 
necessary background on Selberg's method, and its particular application 
to what we call the counting kernel; outline the proof of Theorem [H fill 
in details; and discuss the implications for linear equivalence and strong 
isospectrality. Then we'll given examples to show show how the Theorem 
breaks down in the flat case; show how the Selberg formula works in practice; 
and finish up by discussing some conjectures. 

3 Examples 

In this section, we give examples to show that the trade-offs between bound- 
ary and interior features that are allowed for in the statement of Theorem 
[1] do actually take place: Boundary corners on one orbifold may appear as 
conepoints on the other, while boundary geodesies may migrate to the inte- 
rior. 

All the examples here will be obtained by glueing bunches of congruent 
hyperbolic triangles. The glueing patterns arise from transplantable pairs, 
as described by Buser et al. [2J. In Appendix A below we reproduce a large 
catalog of such transplantable pairs, as computed by John Conway. We will 
refer to this catalog to identify specific glueing patterns. 

Examples of trading boundary and interior geodesies abound. A variation 
on the famous example of Gordon, Webb, and Wolpert [8J yields a pair 
of planar hyperbolic 2-orbifolds of types *224236 and *224623, shown in 
Figure [H (This is presumably the simplest pair of this kind — see Section 
[T2] below.) Each member of the pair is glued together from 7 copies of a 
so-called 346 triangle: a hyperbolic triangle with angles tc/3, 7r/4, ir/G. The 
glueing pattern appears in Appendix A below as pattern 7(3). 

There is no issue here with geodesies passing through the interior of any 
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Figure 1: The pair *224236, *224623. 
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of the triangles that make up these two orbifolds: These can be matched so 
as to preserve length, orientability, and index of imprimitivity. But when it 
comes to geodesies that run along the edges of the triangles, whether along 
the mirror boundary or in the interior of the orbifold, it is necessary to 
balance boundary geodesies on one orbifold against interior geodesies on the 
other, as provided for in Theorem [T] To see this, look at Figure [H and count 
geodesies on the two sides. You'll find that getting the count right is tricky, 
but fun. Beware that boundary geodesies turn back at corners of even order, 
but continue along around the boundary at corners of odd order. Beware 
also of the way interior geodesies bounce when they hit the boundary. The 
answers are indicated in Table [TJ The names 'recto' and 'verso' are short 
for 'orientation-preserving' and 'orientation-reversing', in analogy with the 
names for the front and back of a printed page. The table only shows lengths 
for which there is at least one imprimitive geodesic. Trade-offs between 
boundary and interior geodesies continue at multiples of these lengths. 

Examples of trading corners for conepoints are not as ominpresent as ex- 
amples of trading boundary and interior geodesies, but they are still plentiful. 
Figures |5] and E] show how to construct a pair 6 * 2232233,23 * 22366 from 
transplantable pair ll5'(3) of Appendix A. Note how the order-6 conepoint 
moves to the boundary going one way, while the order-2 and order-3 cone- 
points move to the boundary going the other way. Many other examples of 
cone-trading can be produced using the diagrams of Appendix A, some of 
them much readier to hand than this one. 

4 Background 

Huber [0] proved the result of Theorem [T] in the case of orientable hyperbolic 
2-manifolds. Huber used what would nowadays be seen as a version of the 
Selberg trace formula, which allows us to read off the lengths of geodesies 
from the spectrum in a straight-forward way. Doyle and Rossetti j5] extended 
the result to non-orientable hyperbolic 2-manifolds. In this case we can't 
simply read off the data about geodesies using the trace formula, because 
of interference between the spectral contributions of orientation-preserving 
and orientation-reversing geodesies of the same length. However, it turns out 
that any possible scenario for matching spectral contributions would require 
too many geodesies. 

What we will show here is that, as we indicated in [5j, it is a short step 
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Geodesies for *224236: 



boundary interior total 



length 


recto verso 


recto verso 


recto verso 


2c 


2 2 




■A 3 
2 2 


2a + 26 


1 1 

2 2 


1 1 


3 3 
2 2 


4c 


3 1 3 1 
2*2 2*2 


1 


7 3 
4 4 


4a + 4& 


1,1 1 1 
2~'~2'2 2~'~2'2 


2 • i 

^ 2 


7 3 
4 4 




Geodesies for *224623: 






boundary 


interior 


total 



length 


recto 


verso 


recto verso 


recto 


verso 


2c 


1 

2 


1 

2 


1 1 


3 
2 


3 
2 


2a + 26 


3 
2 


3 
2 




3 
2 


3 
2 


4c 


1 + 1.1 


1 + 1.1 


2 ■ i 


7 


3 


2^2 2 


2^2 2 


^ 2 


4 


4 


4a + 46 


3 1 
2 ' 2 


3 1 
2 ' 2 


1 


7 
4 


3 
4 



Table 1: Gounting geodesies. 




Figure 2: The transplantable pair llg{3). 
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Figure 3: Replace the prototype equilateral triangle in Figure |2] by a 366 
triangle so that six 3-vertices come together in the left-hand diagram. This 
yields a Laplace-isospectral pair of hyperbolic 2-orbifolds of types 6 * 2232233 
and 23 * 22366. 
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from non-orient able surfaces to general orbifolds. The reason is that the 
Selberg formula permits us to read off the data about orbifold features, just 
as Huber read off the data about geodesies in the case of orientable manifolds. 
Then we have only to check that there is still no scenario for matching the 
spectral contributions of the geodesies. 

There are other ways to approach showing that the Laplace spectrum 
determines the data about orbifold features. By looking at the wave trace, 
Dryden and Strohmaier [7] proved the result of Theorem [1] for orientable 
hyperbolic 2-orbifolds. While they did not consider non-orient able surfaces 
or orbifolds, it seems likely that wave techniques could be used to show 
that the Laplace spectrum determines all the orbifold data. The argument 
would be essentially equivalent to the argument we give here, only more 
complicated. 

Another possible approach is via the heat equation. By looking at short- 
time asymptotics of the heat trace, Dryden et al. [6] got information about 
the singular set of general orbifolds (for example, the volume of the reflecting 
boundary). Their results yield information about orbifolds of variable curva- 
ture, and in any dimension. Restricted to hyperbolic 2-orbifolds, the results 
they state don't yield complete information about the singular set. All this 
information is there in the short-time asymptotics of the heat trace, how- 
ever, and presumably it could be extracted using their approach, by looking 
at higher and higher terms in the asymptotic expansion. 

The beauty of the wave and heat approaches is that they can work in great 
generality. The Selberg method depends on having spaces whose underlying 
geometry is homogeneous, such as manifolds and orbifolds of constant cur- 
vature. For studying such spaces, it is a good bet that the Selberg method 
will beat the wave and heat approaches. Of course, the Selberg method can 
be used to treat the heat and wave kernels; the bet is that you will do better 
to consider a simpler kernel, like the counting kernel that we use here. 

5 What we need from Selberg 

Here we assemble what we will need from Selberg for the proof of Theorem 
[U All the ideas come from Selberg |13j . 

Let G = Isom(if2) be the group of isometrics of H^. Note that G has 
two components, corresponding to orientation-preserving and orientation- 
reversing isometrics. A hyperbolic 2-orbifold can be written as a union of 
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quotients of i?^ by discrete cocompact subgroups Fj C G, one for each 
connected component of M: 



We'll denote by F{rj) a fundamental domain for Tj. 

In what follows, we could assume that M is connected, and write 



because the extension to the disconnected case presents no difficulties. We 
prefer not to do this, in part to combat the common prejudice against dis- 
connected spaces. 

Define the counting kernel on H'^ by 



It tells when the hyperbolic distance d{x, y) is at most s. Define the counting 
trace 



The counting trace tells (after dividing by the volume of M) the average 
number of broken geodesic loops on M of length at most s. 

The reason for the name 'counting trace' is that formally, C{s) is the 
trace of the linear operator Lg whose kernel is the counting kernel c(x, y; s) 
pushed down to M. Lg associates to a function f : M ^ R the function 
whose value at x G M is the integral over a ball of radius s of the lift of / to 
the universal cover H^. Lg is not a actually trace class operator, because of 
the discontinuity of the counting kernel, but our expression for C (s) is still 
well-defined. 

We will need the following: 

Proposition 1. The Laplace spectrum determines the counting trace. 




d{x,y) < s 
d{x,y) > s 
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Proof. This is a standard application of the methods of Selberg [T3] . 
which by now are 'classical'. If you aren't familiar with these methods, you 
should look at [13], or failing that, look at section [6] below. The only wrinkle 
here is that the kernel c(a;, s) is not continuous, so the Selberg trace formula 
doesn't apply directly. One way to deal with this is to approximate by a 
sequence of smooth positive kernels that approach the counting kernel from 
above; use the trace formula to see that the Laplace spectrum determines 
the trace of each approximation; and then take a limit, bearing in mind that 
as a function of S, C{s) is positive, increasing, and right-continuous. □ 

The other thing we will need is the 'other side' of the trace formula: which 
will allow us to isolate the contributions to the counting trace arising from 
reflectors and conepoints. 

Let Cl(7, r) denote the conjugacy class of 7 in F, and let Z(7,r) denote 
the centralizer. We have the usual one-to-one correspondence between the 
set of cosets Z(7,r)\r and the conjugacy class Cl(7,r), where to the coset 
Z(7,r)5 we associate the conjugate 6~^'y6. 

Proposition 2. The counting trace can be expressed as a sum of contribu- 
tions from conjugacy classes ofT: 

^(^) = E E Vol({a;GF(Z(7,r,)):rf(x,7a:)<s}). 

i ci(7,r,) 

Proof. It's possible to visualize how this works — for a warm-up, think 
about the case of 1-orbifolds and fiat 2-orbifolds. The trick is to group terms 
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belonging to the same conjugacy class: 




= ^ ^ ^ / c{x,S ^-i5]s)dx 
j ci(7,r,)z(7,r,)5^({,^) 

= ^ ^ ^ / c{5x,jS;s)dx 
j ci(7,r,)z(7,r,)5^(f ,) 

j ci(7,r,)z(7,r,)5^^J.^.) 

= ^ ^ J c{x,'yx;s)dx 
j ci(7,r,V(z(^_r,)) 

= E E Vol({a;GF(Z(7,r,)):rf(x,7a;)<s}). □ 
i ci(7,r,) 

The virtue of this proposition is that we can evaluate the summands using 
simple hyperbolic trigonometry. 

6 The counting trace and the Laplace spec- 
trum 

Here we extract from Selberg [13] the beautiful ideas behind the fact, so 
briefly disposed of in the last section, that the Laplace spectrum determines 
the counting trace. 

We start by going in the opposite direction: 

Proposition 3. The counting trace determines the Laplace spectrum. 

Proof. Denote the heat kernel on by k{x,y;s) = h{d{x,y),t). The 
trace of the heat kernel on M is the Laplace transform of the Laplace spec- 
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trum: 



= E 



Wc can write the heat trace K{t) in terms of the counting trace C(s) by 
means of a Stieltjes integrah 



K{t) = J h{s,t)dC{s). 



This is nothing more than Kelvin's method of images. So the counting trace 
determines the heat trace, and hence the Laplace spectrum. □ 

What we will need here is the converse: The Laplace spectrum determines 
the counting trace. This is hardly surprising — indeed, it would be astonish- 
ing if it did not do so. That would be like discovering that you could not 
determine a well-behaved initial temperature distribution on a semi-infinite 
bar by keeping track of the temperature at the end of the bar. Heat at nearby 
points makes itself felt sooner, and this should be enough to allow us to solve 
the inverse problem, at least in principle. This isn't how our proof is going 
to work, and we don't know if it is possible to prove the result this way. But 
it should be. 

The beauty and power of Sclbcrg's approach is that it gives us a formula 
for how the Laplace spectrum determines the counting trace. But bear in 
mind that we need only the fact, not the formula. 

Here, briefly, is how it works. 

Consider the operator Lg that associates to a function f : M ^ R the 
function whose value at a; G M is the integral over a ball of radius s of the 
lift of / to the universal cover H^. This operator has kernel c{x, y; s), with s 
fixed. 

Ls commutes with the Laplacian, so we can expect to be able to choose 
a basis of eigenfunctions 0o, 0i, . . . of the Laplacian on M that are simulta- 
neously eigenfunction for L^: 

Ls4>i = IJ'i4>i- 
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(If you are worried that the kernel is not sufficiently smooth, approximate 
by a smooth function.) Selberg's great observation is that any eigenfunction 
of A on is automatically an eigenfunction for Ls'. If 

then 

L,(/. = A(A)(/., 

where A is a function that is fixed, known in advance, and independent of M. 
(We'll give a formula for it presently.) This means that any eigenbasis for A 
on M is automatically an eigenbasis for L^, and that the Laplace spectrum 
determines the spectrum (and hence the trace) of Lg. 

Lemma 1 (Selberg). For any A and any Xq G H"^, there is a unique eigen- 
function LOx^xoi^) of A. that is radially symmetric about Xq and normalized so 
that uJx,^o{xo) = 1. 

Proof. Write 

i^X,xo{x) = f(d(x,Xo)), 

and 

/^uJx,xo{x) = Lf{d{x,xo)), 

where L is some second order differential operator that we won't bother to 
write down. Then we have Lf = A/; /(O) = 1; /'(O) = 0. This determines / 
uniquely. (It's some kind of Bessel function.) □ 
Let 

9s{X) = [L,.WA,x-o](a^o), 

and note that this doesn't depend on xq. (It's determined by the function 
we called / in the proof above.) 

Theorem 2 (Selberg). If (f) is any eigenfunction of A on with eigenvalue 
X, then 

Note. While we've formulated this result for operators derived from 
the counting kernel, it applies to any kernel k{x,y) depending only on the 
distance d{x, y). 
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Proof. Fix xq. Denote by the result of making radially symmetric 
about Xq by averaging its rotations. This symmetrization process commutes 
with A, so 

and hence by the Lemma 

But symmetrization about Xq does not change the integral of over a ball 
about Xq, so 

Ls(f){xo) = Ls(t)^^{xo) 

= [Ls(p{xo)ux,xo]{xo) 

= (f){Xo)[LsUJx,xo]iXo) 

= e,{x)(i){xo). □ 

Note. In light of result, we may say that Lg is 'a function of the Lapla- 
cian', and write Ls = Os{A). 

This completes the background for the proof that the Laplace spectrum 
determines the counting trace (Proposition [1]). As we observed there, the 
kernel c(x, y; s) is discontinuous, so the operators Ls aren't trace class, and 
the formal series for C{s) is not guaranteed to be absolutely convergent. To 
get around this, we suggested applying the trace formula to a sequence of 
smooth approximations to c{x,y] s). 

An alternative approach is to integrate the counting kernel c{x, y; s) with 
respect to s to get the continuous kernel 

7/ \ f s — d(x,y) d(x,y) < s 
d(x,y; s) = <^ „ m i~ 

^ ^ \ 0, d{x,y) > s 

The trace of this kernel is 

D{s) = f C{s). 

Js=0 

Since d{x, y; s) is still not smooth, the corresponding trace formula series is 
still not guaranteed to be absolutely convergent. But now we can approxi- 
mate by smoother kernels whose traces will converge pointwise to the value 
of D{s). (This amounts to an alternative way of summing the series for D{s): 
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See Selberg [131 P- 68].) So the Laplace spectrum determines -D(s), and hence 
also C{s). 

Now, in practice the formal sum for C{s) might converge pretty well with 
perhaps a little coddling. At least that is the case for Euclidean orbifolds — 
see Section [11] below. 

7 Outline of the proof 

The counting trace C{s) is built up of contributions from the conjugacy 
classes of Isom(if^). The contribution from the identity is independent of 
s, and equal to the volume of M. The remaining conjugacy classes are of 
four kinds: reflections, rotations, translations, and glide reflections. Their 
contributions all vanish for s = 0, which means that C(0) is the volume of 
M. Translations and glide reflections make no contribution for s smaller than 
the length of the shortest closed geodesic, so for small s the only contributions 
are from the identity, reflections, and rotations. 

Begin by reading off the volume C(0), and subtracting this constant from 
C{s). In what is left, reflections make a contribution of first order in s, 
proportional to the total length of the mirror boundary. Rotations contribute 
only to second order, so we can read off the length of the mirror boundary, 
subtract out its contribution, and what is left is (for small s) entirely the 
result of rotations. 

A simple computation now shows that the contributions of rotations 
through different angles are linearly independent. This allows us to read 
off the orders of the conepoints, and subtract out their contributions. Once 
we've zapped the contributions of the identity, reflections, and rotations, 
what is left of the counting trace vanishes identically for small s, and overall 
is entirely the result of translations and glide reflections, corresponding in 
the quotient to closed geodesies. 

Here we run into the main difficulty in the proof, which is that the con- 
tributions of orientation-preserving and orientation-reversing geodesies of a 
given length are not linearly independent. Indeed, they are proportional, with 
a constant of proportionality depending on the length of the geodesic. So we 
can't simply read off these lengths. Fortunately, we have already established 
that the spectrum determines the lengths and twists of geodesies for hyper- 
bolic 2-manifolds, (See Doyle and Rossetti [5J.) and the argument carries over 
directly to the orbifold case, providing we take care to count geodesies along 
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Figure 4: Possible quotients for the centralizer of a reflection, 
the boundary as half orientation-preserving and half orientation-reversing. 

8 Details 

Denote by R the combined length of all reflecting boundaries. 

Lemma 2. The combined contribution of the reflecting boundaries to the 
counting trace C{s) is 

R sinh - . 
2 

Proof. Let p G Tj be a reflection. There are two possibilities for 
Z{p,rj)\H^: a funnel or a planar domain. (See Figure |H) In either case, the 
contribution to C{s) is i?pSinh|, where Rp is the portion of R attributable 
to Cl(p, Fj), because this measures points in the quotient that are within | 
of the axis of p, and hence within s of their images under the reflection p. 
(See Figure m) Summing over the conjugacy classes of reflections in all the 
groups Fj gives i?sinh |. □ 

Now we turn to conjugacy classes of rotations. Each such is associated to 
a unique conepoint or boundary corner. From a conepoint of order n we get 
n — 1 rotations through angles k = 1, . . . ,n — 1. For each such rotation 
7 the centralizer Z(7,Fj) is cyclic of order n, consisting of just these n — 1 
rotations, together with the identity. The quotient Z(7,Fj) is (surprise!) an 
infinite cone of cone angle — . For a boundary corner of angle - the centralizer 
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Figure 5: The contribution of a reflection to the counting trace. 
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is a dihedral group of order 2n, and the quotient is an infinite sector of angle 
^. This quotient is half of a cone of angle ^, and it should be clear that the 
contribution to C{s) of rotations associated to a boundary corner is just half 
what it is for a conepoint — which is why a boundary corner counts as half 
a cone-point. So we can forget about boundary corners, and consider just 
conepoint s. 

Denote by 5'„(s) the combined contribution to the counting trace of a 
conepoint of order n. This contribution is the sum of the contributions of the 
n — 1 non-trivial conjugacy classes of rotations associated to the conepoint. 
We denote these individual contributions by gk,n{s): 



9n{s) = ^gk,n{s) 

Lemma 3. 



k=l 



gk.n{s) = -/2:rfc(s) 

n " 




where 



Proof. It should be clear that 

gk,n{s) = -f^{s), 

n " 

where fe^s) is the area 27r(coshr — 1) of a hyperbolic disk of such a radius 
r that a chord of length s subtends angle 6. (A fundamental domain for the 
cone hits a fraction ^ of any disk about the center of rotation.) We just have 
to make sure that we have the correct formula for fg. The quantities r,6,s 
satisfy 

sinh I 
sinhr = ^. 



(See Figure El) Thus 

fe{s) = 27r(coshr-l) 



271 (^1 + sinhV - 1^ 



sinh 



2 s 
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Figure 6: The circle whose chord of length s subtends angle 6 
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Hopefully it will seem absolutely incredible that there might be any non- 
trivial linear relation between these functions fg. There isn't: 

Lemma 4. The functions fojO < 9 < tt are linearly independent: No non- 
trivial linear combination of any finite subset of these functions vanishes 
identically on any interval < s < S. 

Proof. It suffices to show linear independence of the family of functions 
\/l + au — 1, 1 < a < oo, on all intervals < u < U. (Set u = sinh^ | 
and a = . \ o .) Actually functions of this form remain independent when 

sm 5 

a is allowed to range throughout the complex plane, and not just over the 
specified real interval. A quick way to see this is to write in Taylor series 

g{u) = VT+^ -l^biu+^u^ + ^ + ... 

and observe (look up?) that only the constant term vanishes: 

61,62,... ^0. 

If a linear combination Y17=i Cig{aiu) is to vanish on < u < [/, then all its 
derivatives must vanish at m = 0: 

n 

^Qaf6fe = 0, A; = 1,2,.... 

j=i 

(We ignore the constant term because all the functions involved here vanish 
at w = 0.) Dividing by b^ gives 

n 

J]c,af = 0, A; = l,2,.... 

i=l 

Could this system of equations for ci, . . . , have a non-trivial solution? If 
so, then the subsystem consisting of only the first n of these equations would 
have a non-trivial solution: 

n 

^Cia*^ = 0, k = l,...,n. 

i=l 

Let di = CiGi, so that 

n 
i=l 
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The matrix of this system of hnear equations is a Vandermonde matrix, 

with determinant nj<j('^« ~ ^j)^ so as long as the are distinct the system 
has no non-trivial solution. □ 

Note. In appealing to the fact that none of the derivatives of \/l + u happen 
to vanish at u = 0, we are seizing upon an accidental feature of the problem. 
We should give a more robust proof. 

Proof of Theorem [T]. Subtract from C{s) from the counting trace the 
volume C(0) and the contributions of the reflectors. Because of the linear 
independence of the functions ftheta, we can identify the number of cone- 
points of highest order. (It should go without saying that this includes the 
half-conepoints at boundary corners.) Subtract out their combined contri- 
bution from C{s). Now we can identify the number of conepoints of the 
next highest order, and subtract out their contribution. Proceed until all 
conepoint contributions have been removed. What remains arises from the 
closed geodesies. In the manifold case, it is shown that the characteristics of 
the closed geodesies are determined by the counting trace, and hence by the 
spectrum. The proof carries over here essentially without change. □ 

9 Implications 

Caveat. The strange and infelicitous notation we use here arises from trying 
to handle disconnected spaces. This leads to things like considering the frame 
bundle of a connected orientable space to be disconnected. There could be 
mistakes in the notation, though hopefully not in the ideas it is meant to 
express. 

Let G = Isom(if^) denote the group of isometrics of H^. Note that G 
has two connected components, corresponding to orientation-preserving and 
orientation-reversing isometrics. A hyperbolic 2-orbifold M can be written 
as a union of quotients of by discrete cocompact subgroups Tj C G, one 
for each connected component of M: 

M = UjTj\H^. 

M is naturally covered by 

M = u, r,\G. 

M is the bundle of all orthonormal frames of M, not just those that have a 
particular orientation; every point x of M is covered by two disjoint circles 
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in M, corresponding to the two orientation classes of frames of the tangent 
space T,j.M . Note that there wiU be two connected components of M for 
each orientable component of M, and one for each non-orientable component, 
because dragging a frame around an orientatation-reversing loop (e.g. one 
that simply bumps off a mirror boundary) will take you from one orientation 
class of frames to the other. 

Note. Even if M happens to be orientable, we should not think of it 
as oriented, because we are discussing manifolds up to isometry, not up to 
orientation-preserving isometry. The Laplace spectrum cannot not detect 
orientation, so it makes no sense to discuss oriented manifolds, and since we 
can't orient our manifolds, there is no obvious good reason to consider only 
orientable manifolds, as has commonly been done in the past when discussing 
isospectrality. But then again, there is no obvious good reason to consider 
only manifolds, as opposed to general orbifolds — possibly disconnected. 

G acts naturally on the right on M, and hence on I?{M). This hnear 
representation of G is analogous to the matrix representation of a finite per- 
mutation group. A finite permutation representation p is determined up to 
linear equivalence by its character x, with x{9) = trp(5f) counting the fixed 
points of g. The situation here should be exactly analogous, the only ques- 
tion being exactly how to define the character. The answer comes from the 
Selberg trace formula. 

For now let's extend the discussion to an arbitrary unimodular Lie group 
G, possibly disconnected. Let F < G be a discrete subgroup with com- 
pact quotient M = r\G. (Notice that we make no mention here of M, 
though in the intended application M will arise from a homogeneous quo- 
tient M — r\G/K.) Denote by 01(^1, G) the conjugacy class of g in G, and 
by Z(g,G) the centralizer of g in G. Introduce Haar measure on Zdg), 
normalized in a consistent (i.e., G-equivariant) way. Attribute to Cl(7, F) the 
weight p^(Z(7, F)\Z(7, G), and define the character Xm to be the function 
associating to g e G the total weight of all conjugacy classes Cl(7, F) for 
which Cl(7, G) = Cl{g, G). Extend the definition of Xm to M = U^- Fj\G by 
linearity. 

As in the finite case, the character is a central function, which means that 
Xuiyh) = XSiihg) for all g,h E G, or what is the same, Xm{9) depends only 
on the conjugacy class G\{g, G). In analogy to the finite case, we can think 
of Xm{9) ^ measuring (in appropriate units) the size of the fixed point set 
{x & M : xg — x}. 
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Proposition 4 (Berard). The character Xm determines the representation 
L'^{M) up to linear equivalence. 

Proof. Using a version of Selberg's trace formula we can write tlie trace of 
the integral operator associated to any smooth function of compact support 
on G in terms of the values of the character Xm- (Cf- Wallach [T51 Theorem 
2.1], Selberg [131 (2-10) on p. 66].) By standard representation theory, these 
traces determine the representation. For details, see Berard [1]. □ 

Proposition 5 (DeTurck and Gordon). The character Xm determines the 
trace of any natural operator on any natural vector bundle over M. 

Proof. Selberg again. (Cf. DeTurck and Gordon PJ.) □ 

With this preparation, we have the following corollaries of Theorem [H 

Corollary 1 (Character is determined). The Laplace spectrum of a hyper- 
bolic 2-orbifold M determines the character Xm- ^ 

Corollary 2 (Representation-equivalence). Laplace-isospectral hyperbolic 2- 
orbifolds determine equivalent representations o/Isom(if^). □ 

Corollary 3 (Strong isospectrality) . // two compact hyperbolic 2-orbifolds 
are Laplace-isospectral then they have the same spectrum for any natural 
operator acting on sections of any natural bundle. □ 

10 Counterexamples 

The analog of Theorem [T] fails for flat 2-orbifolds. In the connected case it 
does go through more or less by accident, just because there aren't many con- 
nected fiat 2-orbifolds. But there are counterexamples among disconnected 
flat 2-orbifolds. We described such examples briefly in |5], and we want to 
recall them briefly here. Someday we will write about them at greater length. 

Here's our premier example, which we call the 1/2 + 1/6 = 2/3 example. 
Let Hi denote the standard hexagonal flat torus ((1, 0), (— |, ^))\R^. Hi 
has as 2-, 3-, and 6- fold quotients a 2222 orbifold H2 (this is a regular tetra- 
hedron); a 333 orbifold H^; and a 236 orbifold H^. (See Figures [7] and [HI) 
Spectrally, 

H2 + He = 2iJ3, (1) 
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r2,r TA^ 1 

\ ATA ATA ATA ATA 
f VAT VAT VAT VaH 

TA ATA ATA ATA A 
f'Af VAT VAT VAT > 
\ ATA ATA ATA ATA 
f VAT VAT VAT VaH 

TA ATA ATA ATA A 
rAT VAT VAT VaT > 

\ /iTa /iTa /iTa /1T7 



r^r y^r T^r i 

V< ATA ATA ATA ATA 
f ViT V^T V^T VA1 

TA ATA ATA AT A A 
ViT ViT ViT > 

k ATA ATA ATA ATA 
f ViT ViT ViT \^ 

TA ATA ATA AT A A 
^AT VAT VAT VAT > 

k iTA iTA iTA iT"^ 



Figure 7: The universal covers of the standard hexagonal torus Hi and its 
quotient orbifolds H^, and Hq. 
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TA 'A AJA 

VaT V; \if V. 

Figure 8: The orbifolds Hi = Q, H2 = 2222, H3 = 333, and Hq = 236. 
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meaning that these spaces have the same Laplace spectrum. (In fact, as 
we'll discuss below, they are isospectral for the Laplacian on fc-forms for 
k = 0,1, 2.) Of course these two spaces match as to volume (1/2+1/6=2/3) 
and number of components (1+1=2). But conepoints do not match. On the 
left we have a 2222 and a 236, so five conepoints of order 2; one of order 3; 
one of order 6. On the right we have two 333s, so nine conepoints of order 3. 

The reason this example is possible is that in the flat case, the contri- 
butions of rotations to the counting kernel differ only by a multiplicative 
constant. Lumping together the contribution of all the rotations associated 
with a conepoint of order n, we get something proportional to (This 
nice simple formula was discovered by Dryden et al. |H].) For a general 2- 
orbifold with variable curvature, it measures the contribution of conepoints to 
the short-time asymptotics of the heat trace. In the fiat case, the short-time 
asymptotics determine the entire contribution, because the contributions of 
all rotations are exactly proportional. 

Combining contributions of all conepoints belonging to connected orbifold 
of various kinds, we get the following totals: 

2222 : 4 ■ I = 6 

333 : 3 • I = 8 

244 : I + 2 ■ f = 9 

236 : I + I + f = 10. 

The isospectrality H2+Hq = 2H-^ arises because 6+10 = 2-8 makes conepoint 
contributions match; we've already observed that the volumes match; and 
geodesic contributions match because on both sides the covering manifolds 
are 2Hi. 

We get a second isospectrality 

Hi + + Hq = 3H2 

because 8 + 10 = 3- 6, 1 + | + | = 3-|, and on both sides the covering 
manifolds are 3Hi. 

Combining these two relations yields other Laplace-isospectral pairs, e.g.: 

Hi + 3Hs = AH2] 

Hi + 4H, = 5H,; 
2Hi + 3iJ6 = 5/^2 • 
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By linearity, all these relations satisfy the condition of having equal vol- 
ume, number of components, and contributions from conepoints. But look 
here: We're talking, in effect, about formal combinations of Hi, H2, H3, Hq, 
so we're in a space of dimension 4. We have three linear conditions, but to 
our surprise, the subspace they determine has dimension 2. Our three linear 
conditions are not independent. If we match volume and number of compo- 
nents, agreement of conepoint contributions follows for free. But why? We 
don't know. 

A similar coincidence happens for square tori. Let Ti denote the standard 
square torus Z^\R^. Ti has as 2-, and 4-fold quotients a 2222 orbifold T2 
and a 244 orbifold T4. We're in a 3-dimensional space, so we expect to be 
out of luck when we impose 3-constraints, but in fact we have the relation 

Ti + 2T4 = 3T2. 

We check equality of volume 1 + | = | and number of components 1 + 2 = 3, 
and then we find that equality of conepoint contributions follows for free: 
2 ■ 9 = 3 ■ 6. Again, why? 

Conjecturally, up to Sunada-equivalence the relations we've given for 
hexagonal and square tori are all there are: See Section [T2] below. 

These relations apply only to the spectrum of the Laplacian on functions, 
not 1-forms. When when we apply the Selberg formula to 1-forms, there is 
still interference between spectral contributions of conepoints, which allows 
our 1/2 + 1/6 = 2/3 example to continue to slip through: It is isospectral 
for 1-forms as well as for functions. But we get additional constraints, which 
taken with those coming from the 0-spectrum force there to be the same 
number of torus components on each side, and this wipes out the other 
examples. The reason we have to have the same number of torus components 
is that this is just half the dimension of the space of 'constant' or 'harmonic' 
1-forms, meaning those belonging to the 0-eigenspace of the Laplacian. And 
since we are expecting just one more linear constraint, this must be it. So, no 
additional mystery here, though when you do the computation it still seems 
a little surprising. 

The isospectrality of the 1/2 + 1/6 = 2/3 example does not hold for all 
natural operators, because the left side admits harmonic differentials that 
the right side does not. 
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0.5 1.0 1.5 2.0 

Figure 9: The counting trace for the *333 orbifold with side length 
computed numerically by summing the contributions of the 11626 eigenval- 
ues with wave number k < 1000. To combat the Gibbs phenomenon (also 
called 'ringing'), we have used the windowing function 1 — (A;/A;max)^ to re- 
duce the contribution of the sum from larger eigenvalues. The upper line is 
the counting trace; the middle is the theoretical contributions from orbifold 
features (reflectors and boundary corners); the lower line is the difference, 
attributable to the identity and other translations. Observe the jump at 0; 
the bend at ^; the jumps at 1 and 2; the jump plus bend at y/S. 

11 Putting the Selberg formula to work 

Here we compute the counting trace for some flat orbifolds, using the Selberg 
formula. 

Figure [9] shows a plot of the counting trace of a *333 orbifold, computed 
approximately from the Laplace spectrum. We can think of *333 as an 
equilateral triangle with Neumann boundary conditions. We have chosen the 
side length to be so that the covering torus is a standard hexagonal torus 
with shortest geodesies of length 1. The shortest geodesic in the quotient 
*333 is the 3-bounce geodesic in the middle, which has length Next 
shortest are the torus geodesies of length 1. Then at length a/3 we get the 
second power of the middle geodesic, which is orientation preserving and 
hence a torus geodesic, plus geodesies parallel to the boundary, which count 
as both orientation-preserving and orientation-reversing. 
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Figure 10: The counting trace of an equilateral triangle with side length 
with Dirichlet boundary conditions. 

Figure [TO] shows another plot of the counting kernel as reconstructed 
from the Laplace spectrum. This time the space is an equilateral triangle 
with Dirichlet boundary conditions. Such a space is a generalized orbifold, 
whose functions change sign under some deck transformation (in this case, 
just those which reverse orientation). Kelvin in effect worked in such a space 
when computing the field of a point charge near an infinite conducting planar 
boundary. The counting kernel counts images with sign ±1. Note the we 
see the same jumps and bends as in the Neumann case above, only now the 
bends are downward rather than upward. 

Now let's try computing the counting trace of a nonorbifold triangle, 
using the same formula as for orbifolds. This can be justified by defining the 
counting kernel to be the Green's function for a variant of the wave equation 
that we may call the counting wave equation. Figures [T2] shows Neumann and 
Dirichlet counting trace plots for an equilateral triangle with angles |, |, ^ 
(side lengths 1,1,0). Because we are using fewer eigenvalues than in the 
graphs above, the bends and jumps are less distinct here. It would be great 
to have more eigenvalues, but these are surprisingly hard to come by! The 
611 Neumann eigenvalues and 591 Dirichlet eigenvalues used here were kindly 
supplied by Alex Barnett. 

In a non-orbifold triangle we expect bends in the counting trace graphs 
corresponding to orbits that are diffracted at the non-orbifold vertices. For 
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Figure 11: The Neumann and Dirichlet counting traces for a triangle with 
angles |, |, ^ and sides 1, 1, 0. 611 Neumann and 591 Dirichlet eigenvalues 
courtesy of Alex Barnett. 
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example, here we see a bend at length 2 sin /racvrS, corresponding to the 
orbit that goes back and forth between the ^ vertex and the midpoint of the 
opposite edge. (This is most visible in the Dirichlet graph, but it is present 
in the Neumann graph as well.) You might enjoy identifying the closed orbits 
that are the source of the other bends and jumps that are visible in these 
graphs. 

Note. Three of the four graphs in Figure [T^ appear to be increasing, just 
as they would be in the case of an orbifold. We don't know if these curves re- 
ally are monotonic, and if so, whether this illustrates a general phenomenon. 
To investigate this, we will want to know how a counting wavefront diffracts 
when it encounters a non-orbifold vertex. We could try some numerical ex- 
periments, if we had eigenfunctions to go along with our list of eigenvalues. 
Alternatively, we could try using the analog for counting waves of Sommer- 
feld's explicit formula for the Green's function of the classical wave equation 
in a wedge. 

12 Conjectures 

As we observed in section [TOl we conjecture, or at least hypothesize, that 
the only spectral relations between fiat orbifolds beyond those coming from 
Sunada's method are those relations we have found for hexagonal and square 
tori; of course these come in all sizes. Now, at this point, Peter doesn't 
remember checking orbifolds with reflectors, though he can't imagine why 
he wouldn't have. But, the hard thing here is not going to be dealing with 
relations between orbifolds. There are two real issues. The first is to see 
that any isospectrality breaks apart into commensurable pieces. The second 
is going to be dealing with the case of tori. Specifically, it seems like for tori, 
the only relations should be those implicit in the theory of Hecke operators. 
The simplest example of a Hecke relation is the isospectrality between 'two 1- 
by-ls and a 2-by-2; two 2-by-ls and a roo-by-roo'. Here 'roo' is my (Peter's) 
proposed abbreviation for 'the square root of two'. 

Moving now to the hyperbolic case. Theorem [T] implies that all hyper- 
bolic triangular orbifolds *ahc (e.g. *237) are determined by their Laplace 
spectrum. It seems likely that all hyperbolic quadrilateral and pentagonal 
orbifolds are also determined. Note that in light of Theorem [H we can treat 
this as a purely geometrical question. This question could most likely be set- 
tled by a computer-aided investigation, along the lines of the computer-aided 
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if 




Figure 14: An isometric pair of type *23424, from glueing pattern 7(2). 

proof that flat 3-tori are determined by their Laplace spectrum. We can ex- 
pect complications in the neighborhood of the isometric pairs of quadrilateral 
and pentagonal orbifolds to which the seven-triangle pairs 7(3) and 7(2) of 
Appendix A degenerate when the triangles are made too symmetrical. (See 
Figures and d) 

This is exactly the kind of issue Rich Schwartz and Pat Hooper faced 
down in the proof of the '100 degree theorem' [11] [12], which states that 
any fiat triangle with no angle larger than 100 degrees has a closed billiard 
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Figure 15: An isometric pair of type *242424, from glueing pattern 7(1). 
Note that this orbifold has 3-fold symmetry, as its Conway symbol suggests 
(but does not imply). 

trajectory. So maybe that investigation would be a better guide here than 
the business about 3-tori. We hope that this problem will turn out to be 
much easier than the billiard problem. 

It seems likely that the hexagons of Figure [1] constitute the simplest 
isospectral pair of hyperbolic hexagonal orbifolds, and possible that all such 
pairs arise from the glueing diagram 7(1). In trying to prove this we can ex- 
pect problems in the neighborhood of the isometric pair we get from glueing 
pattern 7(1). (See Figure [T5l) Another problematic pair arises from pattern 
13a(5). 

A similar computer-assisted approach might show that hyperbolic sur- 
faces of genus 2 and 3 are determined by their Laplace spectra, as suggested 
by Bob Brooks. Here again there is an obvious to expect difficulties, where 
a degenerate Sunada pair yields isometric surfaces. 

One more thing: It seems likely that in any dimension, if two hyper- 
bolic (or spherical) orbifolds are isospectral on fc- forms for all k, then they 
are representation-equivalent. Pesce [TO] showed that that representation- 
equivalence follows from isospectrality for all natural self-adjoint elliptic dif- 
ferential operators on natural bundles. In formulating his result, Pesce simply 
says that representation equivalence follow from isospectrality of all natural 
operators; if you allow non-self-adjoint operators, the problem becomes much 
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easier, because representation-equivalence just amounts to isospectrality for 
the regular action of G on functions on G, and a function on G is a section 
of a natural bundle over r\G. Here we are hoping to make do with only the 
Laplacian on forms. What made Peter think this would be doable is that 
fact that for 2-orbifolds, if you know the spectrum on functions and 1-forms, 
you can simply read off the character Xi using an extension of the 'hard side' 
of Selberg's formula to sections of bundles that are not fiat. This works in 
the 2-dimensional case, where it turns out that if you know the spectrum 
of the Laplacian on functions and on 1-forms, you can just 'read off' the 
character without all the difficulty you face when you try to use functions 
alone. (This is the difficulty that never really arose in this paper, since we 
outsourced the crux of the argument to our earlier paper on nonorientable 
hyperbolic 2- manifolds.) The issue comes down to this: Is an eigenfunction 
of the Laplacian on a. k form automatically an eigenfunction for the averag- 
ing operator where you restrict the form to a ball in the universal cover i7" 
and translate the values back to the center of the ball using the exponential 
map? If this is false, as it may well be when > 2, then the problem could 
be harder than Peter had imagined. 



Appendix: Conway's quilts 

In this section we reproduce John Conway's catalog of transplantable pairs 
coming from small projective groups. Conway generated these pairs using 
his theory of quilts, while Peter watched admiringly. 

The 16 pairs of sizes 7, 13, and 15, which are treelike and thus give 
planar isospectral domains, were reported in Buser et al. [2]. In [2] we also 
presented the first of the size 21 pairs, which yields planar isospectral domains 
for certain triangle shapes. (The second of these size 21 pairs can also yield 
planar domains, so despite what some authors have mistakenly concluded, 
is no special role for the number 17 here.) The rest of these examples have 
been languishing on my webpage. 

There is much to be said about these examples and where they come 
from, but as our purpose here is mainly to bring them to light, we make just 
a few remarks. 

A transplantable pair can be thought of as arising from a pair of finite 
permutation actions of the free group on generators a, 6, c. These actions are 
equivalent as linear representations but (if the pair is to be of any use) not 
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as permutation representations. In the examples at hand, each representing 
permutation is an involution, so these representations factor through the 
quotient F = {a,b,c : = b'^ = = 1), the free product of three copies of 
the group of order 2. 

From any transplantable pair we can get other pairs through the process 
of braiding, which amounts to precomposing the permutation representations 
with automorphisms of F, called L and R: 

L : (a, 6, c) h-t- {aba~^ ,a,c); 

R : (a, b, c) (-)■ (a, c, cbc~^). 

Left-braiding a permutation representation p can be viewed as first conju- 
gating p{b) by p(a), i.e. applying the permutation p{a) to each index in 
the cycle representation of p{b), and then switching p{a) with the new p{b). 
Right-braiding is the same, only with c taking over the role of a. 

Pairs of permutations that are equivalent in this way belong to the same 
quilt. We identify pairs that differ only by permuting a, b, c, or by reversing 
the pair. A quilt has extra structure which we are ignoring here: See Conway 
and Hsu This structure makes it easier to understand and enumerate the 
pairs. But the computer has no trouble churning out pairs belonging to the 
same quilt. 

In the diagrams to follow, the permutations corresponding to a, b, c are 
represented by dotted, dashed, and solid lines. Thin lines connect elements 
that are interchanged, while thick lines delimit fixed points. Thick lines in 
the interior of the diagram separate points that are each fixed, rather than 
interchanged. Sometimes thin lines occur on the boundary of the diagram, 
which means that the boundary must be glued up. We have taken care to lay 
out the diagrams so that there is at most one pair of thin boundary lines of 
each type (dotted, dashed, or solid), so that even without the usual glueing 
arrows there is no ambiguity of how the boundary is to be glued up. 

The label below a pair tells the sequence of braiding operations that gets 
you this pair from the first pair of the quilt. The labels read left-to- right, so 
that {L, L, R, L} means two lefts, a right, then another left. In most but not 
all cases, the quilt has been explored in 'left-first search' order. 

To refer to these pairs, we will write 7(1) for the first pair of quilt 7, 13a(5) 
for the fifth pair of quilt 13a, etc. As this is not very canonical, we reserve 
the right to replace the parenthesized number by a list (possibly empty) of 
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Figure 16: Quilt 7 

Ls and Rs, so that for example 13a(5) would become ISaLLRL, and 7(1) 
would become simply 7. 

For historical reasons, the four quilts of size 11 are called 11/, llg, llh, Hi. 
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{L, L, R} 



Figure 20: Quilt llf 
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{R, R} 



Figure 23: Quilt Hi 
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Figure 24: Quilt 21 
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